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We calculate the finite temperature single hole spectral function and the spin dynamic structure 
factor of spinfuU one-dimensional Tomonaga-Luttinger liquid. Analytical expressions are obtained 
for a number of special cases. We also calculate the single hole spectral function of a spin gapped 
Luther- Emery liquid and obtain exact results at the free fermion point Ks = 1/2. These results 
may be applied to the analysis of angle resolved photoemission and neutron scattering experiments 
on quasi-one-dimensional materials. 



X3 
O 

(N 
> 

oo 

o 
o 
o 



o 



X 



I. INTRODUCTION 



The interacting one-dimensional electron gas (IDEG) 
has intrigued physicists since the pioneering work of 
Tomonagatl. Through the collective effort of many of 
them, and especially due to the advancement of the 
bosonization techniquecTcl, a large number of its prop- 
erties have been discovered. 

The past two decades have witnessed an experimental 
effort to identify and study, using various probes, quasi- 
one-dimensional systems. Among the materials studi 
led one finds organic corapounds such as polyacetylenetl 
and the Bechgaard saltaJ as well as several families of 
inorganic quasi-one-dimensional materialsoH (for exam- 
ple, the transition metal bronzes). The edge states in 
the quantum Hall effect constitute yet another realiza- 
tion of one-dimensional interacting systemaJ. Progress 
made in nanoscale fabrication techniques have made it 
possible to artificially define one-dimensionaUchannels in 
two-dimensional electronic heterostructureal3. Most re- 
cently, experimental evidence has emerged for the ex- 
istence of stripe phases in doped antiferromagnets, such 
as the copper-oxide high temperature superconductorsllj. 
In these phases the doped holes segregate into quasi- 
one-dimensional metallic regions embedded in a predom- 
inantly antifferomagnetic background. 

In recent years angle resolved photoemission spec- 
troscopy (ARPES) has matured into a powerful ex- 
perimental method for probing the single particle 
properties of strongly correlated systems. ARPES 
measurements on the pEjechgaard saltstHi, the transi- 
tion metalchalcogenideal3, the blue bronzesO and the 
cupratescJ have provided evidence in favor of correlated 
one-dimensional physics although full agreement with the 
predicted theoretical picture is still lacking. 

The zero temperature correlation functions of the gap- 
less Tomonaga-Luttinger liquid, and in particular the sin- 
gle hole spectral function which is measured by .ARPES, 
have been calculated pp-viously for the spinlessa as well 
as for the spinfuU caseEElEZl. Extending these results to 
finite temperatures is desirable on the following grounds. 
First, real experiments are carried out at finite temper- 



atures and should be contrasted with theoretical pre- 
dictions relevant for such conditions. Secondly, the 
Tomonaga-Luttinger liquid is a quantum critical system 
which consequently exhibits scaling behavior. It is a rare 
example where not only the scaling exponents but indeed 
the entire scaling functions can be computed explicitly. 
In Section II we obtain the scaling form of the single hole 
spectral function and the spin dynamic structure factor 
of the spinfuU Tomonaga-Luttinger liquid. Closed form 
analytical expressions are derived for a number of special 
cases. These results complement previous numerical j 
nations of the finite temperature spectral functional^ 
and allow for easy investigation of their properties in few 
limits which are of physical interest. Integrals of these 
spectral functions are considered as well. 

When backward scattering of electrons becomes rel- 
evant the IDEG develops a_apin gap and is described 
by the Luther-Emery liquidEHI (Umklapp processes may 
create a charge gap). While the spectrum of the Luther- 
Emery liquid can be readily derived at the special "free 
fermion" point Kg — 1/2, the evaluation of electronic cor- 
relation functions there is non-trivial due to their highly 
non-local form in terms of the pscudo- fcrmions. Several 
authors have made progress in this direction, see Refs. 

and in section III we extend and correct their study 
and obtain exact expressions for the single hole spectral 
function at the free fermion point for temperatures much 
smaller than the spin gap. Some technical details are 
relegated to appendices. 



II. THE TOMONAGA-LUTTINGER LIQUID 

A. The model 

The Tomonaga-Luttinger model embodies the low en- 
ergy and long wavelength physics of the IDEG. It is com- 
posed of two branches of left (77 = —1) and right (77 = -1-1) 
moving massless Dirac fermions constructed around the 
left and right Fermi points of the IDEG. We will consider 
the case where the fermions carry spin 1/2 and denote by 
CT = ±1 their spin polarization. In the absence of back- 
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ward scattering and Umklapp processes the Hamiltonian 
density is given by (?i = 1) 

r;,(T— ±1 

+ 1 92,cPr,{x)p-,,{x) 



E 

ri=±l 



9i,cPnix)Pnix) 
92^s S^(a;)S^^(x) + g^^s S^(a;)S^(x) 



where vp is the noninteracting Fermi velocity and 

Pv = Yi'l,a'^v,'y , 

a 



(1) 

(2) 
(3) 



with T being the Pauh matrices. Despite its appearance, 
the last term in (||) does not break the SU{2) spin sym- 
metry, since S^^(x) = lim^r'^x S^(a::) • S^(a:')/2-l-p^(a;)/8; 
the (72, s term is the only term which breaks this symme- 
try. 

In order to compute the single particle properties of 
the model one uses the bosonization identityOtJ 
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(4) 



which expresses the fermionic fields in terms of self- 
dual fields $,,,cr(x) obeying [^^^^{x),^,,' ^ a'{x')] = 
—iTr6ri^n'^a,tT'Sig'n{x — x'). The Klein factors Frf^a sue re- 
sponsible for reproducing the correct anti-commutation 
relations between different fermionic species and a is a 
short distance cutoff that is taken to zero at the end of 
the calculation. Here and throughout the paper we con- 
sider the limit where the size of the system L is taken 
to infinity and correspondingly ignore terms of the or- 
der 1/L (for a discussion of finite size effects in the finite 
temperature case see Refs. ^5[]2^ ). 

The demonstration of the celebrated charge-spin sep- 
aration in the model is facilitated by expressing <&,),o- in 
terms of the bosonic fields (pc , (ps and their conjugated 
momenta dxOc , dxOg 



(5) 



in terms of which the charge and spin densities are given 
by 



p{x) 
S\x) 



E 
E 

n 



p.,,{x) 



(6) 
(7) 



The decomposition (H) also facilitates the diagonalization 
of the Hamiltonian. It becomes a sum of two independent 
pieces describing noninteracting charge and spin density 
waves which are the elementary excitations of the system 



{dx<Pa) 



(8) 



The velocities of the collective modes are 



(9) 



and the parameters K^, which determine the power-law 
behavior of the correlation functions, read 



i TrVF + g4,a - 92,a 
TTVp + gi,a + 32, a 



(10) 



The effective parameters that enter the Hamiltonian de- 
pend on the specifics of the model for which (|l|) is a low 
energy fixed point. We note, however, that spin-rotation 
invariance dictates g2,s = and consequently Kg — 1. 



B. The space-time correlation functions 



The bosonized expression for the fermion field opera- 
tors (|^) and the fact that the theory reduces to that of 
free bosons make it straightforward to derive expressions 
for the space-time response functions (which we denote 
with a tilde). We will focus on the finite temperature 
single hole Green function 



G<(x,i;T) = (V't_,(x,t)^,,.(0,0)) 



(11) 



We compute this Green function, rather than the more 
usual time ordered or retarded Green functions, because 
it includes only the one-hole states. This fact makes it 
relevant to ARPES whose cross section is directly propor- 
tional to the Fourier transform of . Other Green func- 
tions can be easily obtained from it, since for the model 
d) G>{x,t;T) = (^^,,(0,0)Vi,,(x,t)) = G<{-x,-t-T). 

We will also consider the 2fci? component of the trans- 
verse spin dynamic structure factor, which is measured 
by polarized neutron scattering 



5(x,i;r) = (S^,/(x,t)Sffe^(0,0)) 
+ (S^,/(x,t)S^. (0,0)). 



where 



(12) 



(13) 



and the t are the Pauli matrices. 

The Tomonaga-Luttinger liquid (^ is a quantum crit- 
ical system. It also exhibits spin-charge separation. This 
implies a scaling form for the response functions witi 
separate spin and charge pieces. Specifically one findsc2l 
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X ^'^^ ) ~ f ^ ) (14) 

where we introduced the thermal lengths (ks = 1) 



and the exponents 



(15) 



(16) 



defined so that ja — for noninteracting fermions. Since 
the spin and charge sectors are formally invariant under 
separate Lorentz transformations the functions ga also 
spilt into right and left moving parts 

ga{x, t) = h^^^i {rjx - t) h*^ {ijx + t) , (17) 

where 

h^{x) = [~ismh{x + ia)\~'^ . (18) 
Similarly, for the spin correlation function one findtS@ 
1 



S{x,t-T) 



(2™) 



-2ikFX 



a \ I a ^ 



Xt,c Xt,c 



Xt,s Xt,s 



where 

Ca{x, t) ~ hp^ {x — t) hp^ {x + t) . 
Here we introduced the exponents 



(19) 



(20) 



(21) 



We note that the perpendicular component of the spin 
dynamic structure factor 



S^{x,t;T) ^ {S^^ J {x,t)S^,, (0,0)) 



(22) 



is obtained from the result for the transverse part ([l9| ) 
after multiplying it by an overall factor of 1/2 and using 
the exponents 



0-^ 0-^ 
2 ' 2 ' 



(23) 



instead of ( pT| ) . This holds true also for the spectral func- 
tions calculated in the following subsections. Of coarse, 
in the spin rotation invariant case 5^ = 5/2. 



C. The spectral functions 

It is a general feature of the field-theoretic approach 
to this problem that relatively simple expressions are 
obtained for the space-time dependent correlation func- 
tions. For comparison with experiments, however, we 
are typically interested in the Fourier transform of these 
correlation functions. Conceptually, this is simple, and 
indeed, to evaluate the Fourier transform, we need only 
to perform a two dimensional integral. However, it is not 
generally simple to carry out this calculation analytically. 
Below we consider the cases where this can be done. 

As we noted before, the critical nature of the model 
leads to a scaling form for the spectral functions. A sim- 
plifying feature introduced by the spin-charge separation 
is the ability to express these scaling functions as a con- 
volution of spin and charge parts 

/•OO 

G<(fc,(i;T)= / dxdte'^'^^''+'"'^''"^'^G<{x,t;T) 



a \ ^ I a ^ ^ 



{2t:)^Vc \Xt.cJ \Xt.,, 
X / dqdv gc{q,v) gs{k ~ rq.uj - v) 



(24) 

where we introduce the velocity ratio r = Vs/vc, define 
the scaling variables 



Vsk 
ttT 



UJ 



and 



9a{k,uj)^^h^^^i('^^-^jh,^ 2 



- k 



(25) 



(26) 



h^{k), the Fourier transform of h^{x), is evaluated in Ap- 
pendix A, where some of its properties are listed as well. 
For non-integer values of 7 it is given by 



hy{k) = Re 



(2i)^B ( 



(27) 



where B(x, y) is the Beta function. 

In a similar fashion we obtain for the spin susceptibility 

/"OO 

S{k,[j-T)= / dxdie^t('=+2'=^)"-'"*l5(x,i;T) 



2/3c-l 



2/3s-l 



{2it)^Vc \Xt,c J \Xt,s , 

dqdvCc{q, v) Cs{k — rq,Lj — v) 



(28) 



where 
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(29) 



Henceforth k is measured relative to fc^ and 2k p when 
computing and S, respectively. We also define the 
Fourier transform with respect to rjx. This has the ef- 
fect of making k positive outside the Fermi surface and 
negative inside; thus giving the same expression for the 
spectral function of left and right moving holes. 

Here we would like to note the non-commutativity of 
the limits k,tjj ^ Q and T — )■ in calculating the zero 
temperature dc response of the system. Due to the scal- 
ing form of the spectral functions taking the zero tem- 
perature limit first gives a result which is determined by 
the fc, w — > cx) behavior of the scaling functions, while it 
is the fc, (I; — !■ behavior which is relevant in case k or lu 
is set to zero from the outset. 

Further analytic progress in calculating the spectral 
functions can be achieved in the following cases: 



1. The case Vc = Vs 

Since there is no symmetry between the charge and 
spin sectors, we generally expect that the spin and charge 
velocities are different. This greatly complicates the ex- 
plicit calculation of the Fourier transforms, as the model 
is not "Lorentz invariant" under transformations that in- 
volve both the spin and charge sectors. Such an invari- 
ance is restored if Vc — Vs = v. In this case the correla- 
tion functions have the same form as those of a related 
model of spinless electrons, for which the zero temper- 
ature spectcal functions have been calculated by Luther 
and PescheB. 

The major simplification that follows directly from the 
fact that the charge and spin velocities are equal is that 
the two-dimensional Fourier transform reduces to a prod- 
uct of two one-dimensional transforms - one for the right 
moving piece and one for the left moving piece. Specifi- 
cally one finds, with = 7c + Ts, that 



G<ik,io;T) 



1 



27G 



X h 



I to + k 



io - 



7G + 1 



no 



(30) 

A similar analytic form can be obtained for the spin 
correlation function 



S{k,Co-T) 



1 



2(/3s-l) 



, I Cj + k\ ^ uj — k 



(31) 



2. The spin-rotationally invariant case {Ks = 1) 



We already noted that when the system is invariant 
under spin rotations Ks — 1 (7s =0). At this important 
special point there is no mixing between left and right 
moving spin excitations. As a result the expression for 
the hole spectral function is simplified, due to the ap- 
pearance of the factor ho in the spin part of (^), and 
can be expressed as a single integral 



G<{k,u;T) 



xh. 



^1/2 

" 47r3r V At, 
uj ^ k 



27c 



+ il + r)q 



dq hi ^k — 2rq^ 
uj-k 

— (1 - r)q 

(32) 



When Kc = 1 as well, the integral in (^2|) is trivial and 
we obtain 

G<Ck,^;T) = A^^h, ( ^] h, ^'"^ 



27r2T 11 - rl ^ \ 1-r 



r - 1 



(33) 




where (is ^ I3c + Ps- 



v^k/T 



FIG. 1. MDCs ai LU ^0 (left) and EDCs at = (right), 
for a spin rotationally invariant Tomonaga-Luttinger liquid, 
plotted versus Vsk/T and uj/T respectively, with Vc/va = 3 
and a) 7c = 0, b) 7c — 0.25, and c) 7c = 0.5. 
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The case 7s = 7c = is unique in the sense that there 
is no mixing between left and right moving excitations. 
As a result there are severe kinematic constraints on 
that make it non- vanishing, in the zero temperature limit, 
only in a wedge in fc— plane defined by the lines uj = Vgk 
and a; = Wcfc for A; < 0. While Kg — 1 reflects a symmetry 
of the problem it is unlikely that Kc = 1 {7c = 0) is real- 
ized in any interacting system. However, if the effective 
interactions are not too strong then 7c may be small. For 
instance, for the Hubbard modelEZi, even in the U/t — > oo 
limit, 7c = 1/16. For such systems we expect (p3) to be 
qualitatively correct with the exception of the behavior 
of outside the above mentioned wedge. In contrast 
to the result for 7^ = the T = support of G^, for 
non-zero 7c, extends up to the line lu — —vjt with fc > 0. 
Nevertheless, if the mixing between left and right moving 
charge excitations (i.e 7c) is small the amplitude of G^ 
outside the wedge is small and its gross features resemble 
the 7c = result. 

The customary way to present ARPES data is by plot- 
ting momentum distribution curves (MDCs) and energy 
distribution curves (EDCs). These curves are cuts in 
G^{k,uj) of constant lo and constant k respectively. In 
Fig. 1 we present MDCs at the Fermi energy {cu = 0) and 
EDCs at the Fermi wave-vector (fc = 0) for a spin rota- 
tionally invariant Tomonaga-Luttinger liquid for various 
values of the parameter 7c. 



G< (fc, ^; T) = ^ (3^) " h,^^^. i^-k)h. (fc) . 

(36) 

D. Integrals of the spectral functions 

Integrals of the spectral functions, such as the density 
of states and the momentum occupation number, are of 
interest too. These simpler quantities give us a qualita- 
tive view of the spectrum, without too many complicated 
details. 



1. p^iu)) and So{ui) 

The calculation of the density of states p'^{co) and its 
analogous quantity for the spin susceptibilities Sq{lo) is 
straightforward. The results are 

p< {Co; T) = / ^e-*G< (0, t) = / dk G< (fc, ~cj) 

(37) 



3. The case Vs/vc 

In the limit where one of the velocities is much smaller 
than the other the calculation becomes again more 
tractable. In many physical systems Vc > Vg and we 
will consider this case (for example for the t — J model 
away from half filling we have Vg/vc ^ J/t < 1). In the 
limit r = Vg/vc ~* the spin piece in Eq. ( p^ ) becomes 
q-independent and any singularities in G^ disperse with 
the slow velocity Vg ■ The q integral is then readily evalu- 
ated by expressing the factors in gc as Fourier transforms 
of their real-space counterparts. The result is 



G<(fc,cI;;T) 



a \ ^ I a ^ ^ 



/OO 
du /^27c+^ {uj — k — 2u) h^^^i (k + y) h^^ (i^) . 
-00 

(34) 

A similar calculation gives for the spin spectral function 



S{k,Co;T) 



1 



2/3c-l 



2f3s-l 



X / dv h2f3^{u) — k — 2^) hp^{k + ly) hi3^{h') . (35) 



At the spin-rotationally invariant point {Kg — 1) the 
integral in (|3^) can be performed to give a result for G^ 
in the limit Vg/vc — > but arbitrary 7c 



So{^;T) 



^e-*5(0,f) 

2/3,- 



dkS{k, —uj) 



STT-^a ^JVcVg \ At, 



At,; 



2/3. 



h2l3s {-^) ■ 
(38) 



2. n{k) and S{k) for the case Vc = Vs 

The evaluation of the momentum occupation number 
n{k) and the spin structure factor S{k) is complicated 
by the fact that they can not be expressed in terms of 
the functions hy{k). Contrary to the scaling functions 
considered above the scaling functions of these quanti- 
ties depend on the cut-off a~^. It is, however, possible 
to compute them in the case where the charge and spin 
velocities are equal. 



n(fc;T)= / da;e'''('=+'='=-)^G<(x,0) 



diu 

TT \XtJ 



G<(fc,w) 



X Im < / dx 



[cosh(2a;) - cosh(2ia/AT)]T<=^+M ' 

(39) 
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where again 70 = 7c + 7s- Since the integrand in ( |39| ) 
is imaginary in the interval [ia/\T, 0] of the imaginary 
axis, we can add the i nteg ral along this interval to the 
one already present in ( p9|) without affecting n{k). 



n{k]T) 



27G-1 



27G 



"^^"^ 4a/A? [cosh(a;) - cosh(2ia/AT)]^«+i 

(40) 

Formally, for 7^ < 0, Eq. ( ^0| ) is related to the inte- 
gral representation of the associated Legendre function 
(defined with a branch cut from —00 to 1) (3(j(z), see 
Ref. |28[ Since n{k) is regular in 70 we can use analytic 
continuation to obtain for all non-integer 



21G-1/2 / Q 



27G 



X Re<| e-''^T'^[sinh(2m/AT)]~^''''+^/^^ 



(41) 



It can be shown that for 7^ < 1/2 and A; ^ (Ara) 
this expression reduces, in the limit a — > 0, to 

n{k- T) = I^I^Re|r (7G + 1/2) e'^^^^+i/^) 

^\ ^^-^ r(l/2^7G)r[(l + 27G-»fc)/2] 1 
^tJ r[(l - 27G - ifc)/2] ]■ 

(42) 

Using similar manipulations we obtain, for non-integer 
I^S = (ic + (3s, the spin structure factor 



S{kT) 



9/35-3/2 / \ 2/3S-1 



2^ 



S'(fc,u;) 



X Rel e"(i/2-&)[sini,(2ia/AT)]i/2-'^^ 



><Q2^;/-i/2[-«(2«/At)]|, (43) 

which for a -> 0, /S^ < 3/2 and fc < {Xto)-'^/'^ tends to 

^(^; = £lL_MRe J r (/35 - 1/2) e-(i/2-fe) 
47r''/^a 



a_>l''^^-^r(l/2-/35)r(/35-^fc/2)' 
At, 



r(l -I3s- ik/2) 



(44) 



III. THE LUTHER-EMERY LIQUID 

A. The model 

In the Tomonaga-Luttinger liquid both the charge and 
spin excitations are gapless. Spin and charge gaps may 
open up as a result of adding to the Tomonaga-Luttinger 
Hamiltonian (^ terms describing backward and Umk- 
lapp scattering respectively. In the following we will 
assume that the IDEG is sufficiently incommensurate 
so that Umklapp scattering may be neglected. Conse- 
quently the charge sector remains gapless and contin- 
ues to be described by a free bosonic theory [the charge 
part of Eq. (^]. Including a backward scattering term 

51 Sij V'J iV'-,, -i'0r;,-iV'-?),i results, after bosonization, 
in a spin Hamiltonian density of the sine-Gordon type 



Ks 



{2na)' 



■ cos(v87r(?!)s) 



(45) 



The gi perturbation is relevant for Ks < 1, in which 
case a spin gap is-dynamically generated according to the 
scaling relationESi Ag ~ (i!s/a)[(7i/27r^Ws]^/(^~^^=' and 
the excitations are massive spin solitons. In the spin 
gapped phase the problem is most simply treated in terms 
of spin fermion fields and their mode decomposition 
C'± {0) (obeying the usual fermionic anti-commutation re- 
lations) in momentum space 



= F,, exphivW2(^. - 2770s)] 



47rw, 



(46) 
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where we introduce the spin correlation length = 
Vg/As and use the rapidity representation k = sinh^/^^. 

For Kg — 1/2, which is known as the free fermion or 
Luther-Emery pointE3, the refermionized Hamiltonian is 
non-interacting and massive with a gap Ag = gi/2iTa 

))=±i 



deEs{e) \cl{e)c+{e) + cl{e)c^{e) - 1 



(47) 



where the spin excitation spectrum is 



Es{B) = A, cosh0 = + (vskf = Es{k) . (48) 

In the following we will concentrate on computing the 
correlation functions at the Luther- Emery point. We will 
comment briefly on the effects of deviations from this 
point. 



B. The spin part of the spectral functions 

In contrast to the Tomonaga-Luttinger model the cal- 
culation of the spin part of most spectral functions for 
the Luther-Emery liquid is no longer trivial. The diffi- 
culty lies in the fact that generically the refermionized 
form for the spectral functions involves highly non-local 
operators. We start by evaluating the spin contribution 
to the transverse spin form factor, which fortunately has 
a simple representation in terms of the pseudo-fermions. 
We then consider the single hole spectral function that 
belongs to the wider class of functions which do not ad- 
mit such a simple form. Interestingly, it is still possible 
to obtain an exact expression for this function too at the 
Luther-Emery point. This, exact result extjcnds and cor- 
rects earlier work of VoilO and Wiegmannc3. 

In computing the various spectral properties of the sys- 
tem we can distinguish between two temperature regions. 
At temperatures large compared to As the spin gap can 
be ignored, and the results for the Tomonaga-Luttinger 
liquid apply. If the temperature is small compared to 
the spin gap we can evaluate the spin contribution to the 
correlation functions in the zero temperature limit while 
introducing errors of order exp(— A^/T). Henceforth we 
will be concerned with temperatures in the second region 
only. 



1. The spin part of S{k,uj) 

The spin piece of the transverse spin correlation func- 
tion has the following simple form in terms of the spin 
fermion fields 

S,ix,t) = (*I(a:,i)«'Li(a:,i)*-i(0,0)*i(0,0)) . (49) 



Since the theory reduces, at the Luther-Emery point, to a 
theory of free massive fermions, the corresponding spec- 
tral function can be readily computed with the result, for 

T = 0, 



4EUk/2) 



Avl\qiE,{q2) - q2E4qi)\ 



e[~UJ-2Es{k/2)] , 
(50) 



where the spin excitation spectrum Egik) is given by 
Eq.(^8[) and qi 2 are the solutions to the quadratic equa- 
tion LJ + Es{q) + Es{k — q) ~ 0, that is 



_k to 4A2 



(51) 



2. The spin part of {k,uj) 

The refermionized form of the single hole Green func- 
tion (below we consider the case rj = — \ and (7=1) 

Gs{x,t) = (C/I(x,t)^'Li(a:,i)*-i(0,0)C/i(0,0)) , (52) 

4 4 

involves the non-local vertex operators 



Uaix) — exp[iv87ra0s(a;)] , 



with 



(j)s{x) 



(53) 



(54) 



it is evident that and 



From Eqs. (|6D, (^ and (| 
Ua create and destroy a single spin soliton or an integer 
number of soliton pairs, respectively. Therefor, the Green 
function consists of a coherent oae,spin soliton piece and 
an incoherent multi-soliton piecew 

Gsik,Lu) = Z,ik)d[Lo + Esik)] + Gi™"'")(fc,c^) , (55) 

where the multi-soliton piece is proportional (at T = 0) 
to Q[-~Lu — 3Es{k/3)]. Deviations from the Luther-Emery 
point in the case Kg < 1/2 will result in the formation 
of a spin soliton-antisoliton bound state, a "breather" , 
which can shift the threshold energy for multi-soliton ex- 
citations somewhat. 

At the Luther-Emery point the form factors of the ver- 
tex operators, i.e. their matrix elements between the 
vacuum and various multi-soliton states, are known ex- 
actly. This fact enables us to obtain exact results for the 
different parts of the spectral function. The details of the 
calculation are presented in Appendix B. For the spectral 
weight of the coherent piece we find 



Zs{k) 



8c2 

TT 



26 
a 



Wcfc 



Esik) 



(56) 
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where c = 0.101. We|-s5jould like to note here that a 
simple scaling argumentE3 allows us to obtain the depen- 
dence of Zs on for arbitrary Kg < 1. It follows from 
the observation that the sine-Gordon theory is asymptot- 
ically free and hence the dependence of Gs on the short 
distance cutoff a is unaffected by the opening of a spin 
gap. Since in the absence of a gap Gs is proportional to 
(j27,-i/2 jg ^ jjiatter of dimensional analysis to see that 

Zs(k) - (6/a)*-'^Vs(fc6) , (57) 

where /s is an undetermined scaling function. 

The incoherent piece of Gf{k,uj) consists of contribu- 
tions from processes involving intermediate states con- 
taining N = 2n + 1 spin solitons, n = 1,2,3, 
Each such contribution starts at an energy threshold 
-NEsik/N) 



N=3,5,- 



-uj - NE. 



(58) 



At the vicinity of the threshold i.e. for fc^s ^ 1 and 
\uj + NEs{k/N)\/As < 1 the behavior of is given 

by [uj + NEs{k/N)]"^+-^-'^. In particular the 3 soliton 
part reads then 



Gf'°'Hk,Lo) 



_8c^ 
V37r2 



26 



3Es{k/3) 



A? 



(59) 



As discussed above, for large \uj\ the incoherent piece 
should asymptotically approach the Tomonaga-Luttinger 
result w2t=-3/2. 



C. The single hole spectral function 

Once the spin part is calculated there still remains the 
task of convolving it with the corresponding charge part 
in order to obtain an expression for the full spectral func- 
tion. Analytically this is difficult, and like in the case 
of the Tomonaga-Luttinger liquid, progress can be made 
only in special cases. Below we carry out the convolution 
for the single hole spectral function. We consider tem- 
peratures well below the spin gap scale and correspond- 
ingly use the above derived zero temperature results for 
the spin part. For the gapless charge degrees of free- 
dom we continue to use the finite temperature Tomonaga- 
Luttinger expressions. 

The spectral function can be written as a sum of two 
contributions = G1+G2 coming from the convolution 
of the charge part with the coherent (single soliton) and 
incoherent (multi-soliton) pieces of Gf. The frequency 
integral in Gi is readily evaluated with the result 



Giik,Lo;T) 



A2 



-) 



270+1 



26 
a 



dq 



1 



Vs{k - q) 



Esik-q) 
UJ + Es{k - q) - Vcq 



Es{k - q) + Vcq 



2ttT 



2ttT 



(60) 



where here At — Vc/i^T. 

In Fig. 2 we present representative contributions of the 
single spin soliton piece, Gi, to the MDCs and EDCs of 
a Luther-Emery liquid with various values of the charge 
exponent 7c. Using some of the results derived below, we 
also indicate the asymptotic behavior of the three spin 
soliton contribution, G^2 to the EDCs in the vicinity 
of its zero temperature threshold lo = — 3As. 



a j 




b ^ 




c A 








FIG. 2. The single spin soliton contribution, Gi, to the 
MDCs at a; = (left) and EDCs at fc = (right), of a 
Luther-Emery liquid [K^ = 1/2), with Vcjvs = 3, As/T = 3 
and a) 7c = 0, b) 7c = 0.2, and c) 7c = 0.4. The asymp- 



(3 sol) 



to 



totic contribution of the three spin soliton piece, G 
the EDCs, near its zero temperature threshold uj = — 3As, 
is indicated by the dashed lines. This contribution was cal- 
culated using Eq. (^^ in the case 7c — 0. For the EDCs 
with 7c = 0.2, 0.4 we used the asymptotic result for the case 
Vs/vc — > 0, Eq. (Ba). Other multi-soliton processes start to 
contribute near loTAs — 5,7, ■ ■ ■, as discussed in the text. 



We now restrict the discussion to the following special 
circumstances: 



8 



1. The case Vs/vc 







2. the case 7c = 



In this case the remaining integral in ( |60| ) is straightfor- 
ward since the q dependence of the spin part and Es{k—q) 
disappears and we are left with 



4c^ a 



a 



1 - 



^ At 

Es(k) 



26 



uj + Es{k) 



ttT 



(61) 



which in the limit of zero temperature reduces to 



TT r(27, + l/2) 



27c 



26 



[-UJ - E,{k)p^-"^ Q[-uj - E,{k)] . (62) 



For the 3-soliton contribution to G2 we obtain, assum- 
ing fc6 < 1 and \lo + iEs{k/'i) \ < A^, 



G. 



(3 sol) 



{k,oj-T) 



Ac' 



1 



V3^4 A,T VAt 



27c 



26 



uj+?,E,(k/?,) 



v-uj- 3Es{k/3) 



(63) 



which at r = equals 

fi^2 



G 



(3 sol) 



(A:,c^;0) 



V3^2r(27e + 5/2) Vw, 



27c -Hi 



26 



X — [-c^ - 3^,(fc/3)]'^=+^ ehc^ - 3Es{k/3)] . 

(64) 

For this case the N = 2n -I- 1 soliton contribution to 
G2 is proportional, in the vicinity of its threshold, to 
[-UJ - 7Vi;^(fc/iV)]2T=+"'+"-i/2e[-u; - NEsik/N)]. 



It is also possible to derive closed expressions for the 
spectral function when 7^, = but r = v^/vc is arbitrary. 



Gi(fc,w;T) 



4c^ 1 /26 

TT^ 1 — r Vr \ a 



X h 



u! — rvgk + ^ [ruj — VgkY -f (1 — f'^)A^ 



7r(l -r2)r 

ruj — Vgk -\- \J (ruj — Vgk)'' + (1 — r^) A^ 
^{rw - VskY + (1 - r2)Af 



(65) 



Using the asymptotic form (A9) of /i^(fc) one can obtain 
the zero temperature limit of Gi. When r < 1 we find 



Gi(fc,w;0) 



8c2 



r3/2 V 1 



1 + r 



26 
a 



rtj — Wsfc + -^Z [ruj — Vgk)' + (1 — r2)A^ 
y^iruj - i;,fc)2 -f (1 - r2)Af 
e[-a;-£;,(fc)] 



-a; -|- rwsfc — (ruj — Vgk)' + (1 — r^)A^ 

(66) 



The kinematics in the case r > 1 is more involved. In 
particular we find that for sufficiently large r it is possible 
to distribute the total momentum between the spin soli- 
ton and the gapless charge modes in a way that gives a 
contribution at frequencies smaller than Es{k). (In this 
case the effect exists only for fc < 0. This is a special 
property of the point 7c = which precludes mixing of 
left and right moving charge excitations.) 



J 



Gi(fc,a;;0) 



8c2 



71-3/2 y 7. _ 1 

E 



26 



Vgk — ruj — [ruj — Vsk)'' + (1 ^ r2)A2 



yjiruj - VskY + (1 - r2)A2y^Lj - rVgk + (J^J{ruj - VskY + (1 - r2)A2 
cr = 1 if < -Es{k) 

a = -1, 1 if fc6 < and - E,{k) < uj < -^k-Vr^A^ ■ 

I 



For fc6 < 1 and \uj + 3Es{k/3)\ < A^ the 3-soliton 
contribution to G2 reads 



(67) 



VcQ-uj - 3Es[{k - q)/3] 



X e{vcq-Lo-3Es[ik-q)/3]} . (68) 



G'i'°''>{k,iu;T) 



4c2 7^ /26 
V37r3 As \ a 



dq hi {Xti) 



When r < 1 the zero temperature limit of the above is 
readily evaluated giving 
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G 



(3 sol) 



33/27,5/2 
uj + 3Es{k/3) 



a 



Q[-uj -3Esik/3)] . (69) 



The zero temperature limit when r > 1 is, once again, 
more comphcated. However, it is possible to show that 
in the range of vahdity of Eq. (|68|), i.e. for k^s ^ 1, and 
for 1 < r < 3 it coincides with the expression for the case 
r < 1 Eq. (^9|) . Deviations from this behavior may occur 
only if the velocity ratio is large and satisfies r ^ 3. 
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APPENDIX A: THE FUNCTION h-y{k) 

The function h^{k) is real 

/•OO 

h^{k) ^ / dx e''^^ h.y{x) 

dx , . . S> . (Al) 



lim 2Re 



-i sinh(a; + ia)\^ 



Although the imaginary part of the integral diverges as 
*" for 7 > 1, its real part, and hence also h^{k), are 
analytic for all values of 7. Substituting y — e~^^ we find 



h^{k) = lim Re <^ {2iy / dyy 



--i(l-e-2«^^^^ 



(A2) 



The integral is analytic in the limit a — > for 7 < 1. In 
this range the exponential factor in (|A2| ) can be dropped 
and it red,uces to the integral representation of the beta 
functionEl. We can then use analytical continuation to 
obtain for all non-integer values of 7 



h^{k) — Re 



(A3) 



For the integers, hn{k) can be calculated as follows. 
First, we obtain using standard residue technique 



h,{k) = 27r/+(^fc) , 

h2ik) = 2TTkf_{-Kk) , 



(A4) 
(A5) 



where f±{k) = {e^ ± l) ^ are the fermionic and bosonic 
occupation functions. We then integrate (Al) by parts 
twice to find the recursion relation 



l-n+2 



(k) 



which implies 



n{n + 1) 



n-l 



h„{k) 



(A6) 



2tt 



h2n{k) 



27r 
r(2n) 



m=0 

n-l 



k \[ [{2mf + e]f^{nk) . (A7) 



m— 1 



We also note that 



ho{k) 27r(5(fc) 



(A8) 



Finally, the asymptotic behavior of h^{k) for large |fc| 
is easily evaluated with the result 



h,{\k\ 



2ti 



Q{-k){-ky-^ 



(A9) 



Since the spectral functions have a scaling form in the 
variables lu/T and Vak/T Eq. (A£) also determines the 
zero temperature limit of these functions. 



APPENDIX B: CALCULATING Gf{k,uj) AT THE 
LUTHER-EMERY POINT 

Inserting the resolution of the identity into Eq. ( |5^ ) 
one obtains 



Gfix,t) = y2{0\U\{x,t)\N)(N\^l,{x,t)^^,{0,0)\N') 

x(7V'|C/i(0,0)|0) . (Bl) 



N,N' 



The matrix elements of the vertex operators appearing 
in (Bl) are known as the form factors of these operators. 
They have been derived, at the Luther-Emery point (see 
also Ref. |3^), by a variety of Wfys. They were first 
obtained by Schroer and Truong£3 who normal ordered 
the vertpj operator with respect to the spin fermions. 
SmirnovE^ derived them for a = 1/2 using bootstrap ax- 
ioms. Most recently they were calculated-, using mon- 
odromy relations by Bernard and LeClairE3. Since the 
fields J7q.(0, 0) are neutral with respect to the topologi- 
cal U{1) charge of the solitons (±) the form factors are 
non- vanishing only for U{1) neutral states 



|t/„(x,t)C|(02„) • • • Cl{9n+l)CU0n) ■ ■ ■ CUOlM = 



(-1) 



n(Ti-l)/2 



sm na 

2TTi 
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X exp 



■ 2n 

i I — sinh(?fc — A,,t cosh( 

. k=l 



n 
k=l 



<k<j<n 



sinh 



(^)sinh(^ 



ni<fe,,<„cosh(^^) 



(B2) 



where Va, the vacuupi expectation value of the vertex 
operators, is given bji23 



F„EE (0|C/„ (0,0) |0> =c(a) 



26 



with 



c(a) = exp 



dt 
T 



sinh^ (at) 
sinh^ t 



(B3) 



(B4) 



Furthermore, since creates and destroys a single 
soliton the state |A^') can differ from |A^) by 0,±2 soli- 
tons only. It is also easy to check that Gf{k,uj) is real. 
Using Eqs. (^6|) and ( |b^ ) one obtains that each term 
in the Fourier transform of Eq. (Bl) is proportional 
to )/2 times a real expression. Thus only terms 

with {N + N')/2 an even integer should be considered. 
Combining these two observations we conclude that the 
contribution to Gf{k,oj) comes solely from the terms 
N = N' (the amplitudes for the cases N = N' zL 2 are 
finite but cancel each other). 

The coherent piece of the spectral function is due to 
the terms iV = and N = 2. The evaluation of the first 
is straightforward with the result 



1 - 



S[u; + Esik)] , 



where V = Vi = c(i)(26/a)-i/i6. The cutoff depen- 
dence of the above result is in conflict with the general 



scaling argument that was given following Eq. (p6|). It 
represents corrections to scaling coming from irrelevant 
operators. The leading scaling behavior is recovered by 
considering the contribution of the N = 2 term to the 
coherent piece 



1/2 



32713^ 



dp 



cosh - 



d9 e~^^*['^/'^=^ sinh e- As t cosh 6 



The p integral is divergent and a lower cutoff — ln(26/a) 
(corresponding to a cutoff k = — 1/a) should be intro- 
duced. We then find for the leading contribution to the 
coherent piece 



Gi™'"='^)(a;,i) 



Ac' 
xKi 



26 
a 



Vst 



^a;2 - [vst] 



y/x'^ - [vsty + iet 

-I- ie-t 

, (B5) 



where Ki (x) is a modified Bessel function. Fourier trans- 
forming it one obtains Zs{k)6jui+Es{k)] with the spectral 
weight Zs{k) given by Eq. (|5q). 

The 2n -|- 1 soliton contribution to the incoherent piece 
of the spectral function comes from the terms N = 2n 
and N = 2n + 2. As is the case with the coherent piece 
the former contains corrections to scaling while the latter 
is responsible for the leading behavior which is propor- 
tional to 



V' d9i--- d92n+l 



dPe 



T\.n+2<j<2n+l 

sinh 
ni<j<„+i cosh 



2ri+l 



2n+l 



'-,«+! , 2n+l \ n ".inh^ ( Sk-Oj \ TT ■ ,2 / flfc-flj 

[ 1 1 fl 1 •lil<fe<i<"+l V 2 J i.ln+2<k<j<2n+l^^^^ \ 2 



2 ^ 

j=n+2 



n 



n+2<fe<2ri+l , l<j<n+l 



cosh" 



(^) 



The exact evaluation of the above integrals is difficult but 
their behavior near the threshold (2n + l)Es[k/ {2n + 1)] 
may be extracted in the following way. Writing uj — 
-{2n+\)E,\kl(2n^\)\ - Aw one finds that for A:^ < 1 
and Aw ^ A^ the second (5-function restricts the inte- 
gration region over Qi to a small ball near the origin. 
It is then legitimate to expand the integrand to lowest 
order in these variables. The integral over gives (to- 
gether with the factor V'^^ the overall cutoff dependence 
(26/a)'^^'^ in accord with the general scaling argument. 



By changing variables to spherical coordinates it is easy 
to check that the remaining integral over Qi is propor- 
tional to Aw" "'""'~-'^0(Aa;). A detailed evaluation of the 
3 solitons case results in Eq. (p9|). 
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